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Abstract 



The issue of the number, nature and sequence of phase transitions 
C^ ' in the fully frustrated XY (FFXY) model is a highly non trivial one 

due to the complex interplay between its continuous and discrete de- 
grees of freedom. In this contribution we attack such a problem by 
means of a twisted conformal field theory (CFT) approach [1] and 
show how it gives rise to the U (1) Z2 symmetry and to the whole 
spectrum of excitations of the FFXY model [2]. 



1 Introduction: the state of art 

The phase diagram of the FFXY model has been the subject of intensive 
studies in the last thirty years, due to the presence of the mixed symmetry 
U (1) ® Z2. But a full and definitive answer to such a question still lacks; in 
this paper we address the problem by means of a twisted CFT approach [1]. 
The XY model on a square lattice in the presence of an external magnetic 
field transversal to the lattice plane is described by the action: 

^ = —j^^^^^(Vi-Vj-Aj), (1) 

where {ip} are the phase variables on the sites, the sum is over nearest neigh- 
bors, J > is the coupling constant and the bond variables Aij = j^ f^ A-dl 
satisfy the full frustration condition J^piaquette ^ij ~ ^- Choosing the Lan- 
dau gauge in such a way to get a lattice where each plaquette displays 
one antiferromagnetic and three ferromagnetic bonds, we obtain two ground 
states with opposite chiralities. The discrete Z2 symmetry of the FFXY 
model is broken at low temperature and will be restored beyond a certain 
temperature after the formation of domain walls separating islands of op- 
posite chirality. The Ising transition overlaps to a vortex-unbinding tran- 
sition, which is associated with the continuous U{1) symmetry [3]. The 
action (1) can be rewritten as a fractionally charged Coulomb gas (CG) de- 
fined on the dual lattice, Hcg = ~4f '^ry ("^ (^) + 2) ^ (^' ^') ("^ (^') + 2) ' 
where hmi^-r'^oo C (r, r') = log|r — r'| + ^tt and the neutrality condition 
X^r ("^ (^) + i) ~ holds. Such a model exhibits two possible phase transi- 
tions, an Ising and a vortex-unbinding one. The issue whether there are two 
distinct phase transitions, Ty > T^w or Ty < T^w with Ty and T^w mark- 
ing respectively the breaking of U{1) and of Z2 symmetry [4], or a single 
transition with the simultaneous breaking of both symmetries [5] has been 
widely investigated. The model allows for the existence of two topological 
excitations: vortices and domain walls. Vortices are point-like defects such 
that the phase rotates by ±27r in going around them [3]. A domain wall can 
be viewed as a line on the square lattice, each segment of which separates 
two cells with the same chirality [6] . If the domain walls form a right angle, 
such corners must behave as fractional vortices with topological charge ±1/4 
[3]. There exists a temperature T^w such that, when T > Tdw, dissociation of 



bound pairs of fractional vortices is allowed, which triggers the dissociation 
of pairs of ordinary vortices. The system undergoes two phase transitions 
with temperatures such that Ty < T^w [7] . At finite temperatures kinks may 
appear on the domain wall. Simple kinks must behave as fractional vortices 
with topological charge ±1/2 while a double kink does not introduce mis- 
matches in the phase distribution. At low temperatures, all simple kinks are 
bound into neutral pairs. As the temperature increases, a phase transition 
in the gas of logarithmically interacting kinks leads to pair dissociation and 
emergence of free simple kinks [6]. That takes place at Tk < Ty [8] and 
produces two distinct bulk transitions with Ty < Tdw A more general con- 
clusion is reached by studying the coupled XY-Ising model [9, 10], which is in 
the same universality class of the FFXY model. Such a model can be intro- 
duced starting with a system of two XY models coupled through a symmetry 
breaking term [9]: 



H = A 



Y^ Y^ cosfv^«W-V^^'^( 
i=i,2 i^^y^ 



+ /iVcos2(v9('n^)-V^^'H^))- 

(2) 



The limit /i ^ corresponds to a full decoupling of the fields (/9*^*\ i = 
1,2, while the h ^ oo limit corresponds to the phase locking (f^^\r) — 
(p^'^\r) = Tij, j = 1, 2; as a consequence the model gains a symmetry U{1)0Z2 
and its Hamiltonian renormalizes towards the XY-Ising model [9]. Its phase 
diagram [9, 10, 12] is built up with three branches which meet at a multi- 
critical point P. Two branches describe separate Kosterlitz-Thouless (KT) 
and Ising transitions while the third (PT) corresponds to single transitions 
with simultaneous loss of XY and Ising order. It becomes a first order one at 
a tricritical point T and seems to be non-universal [9]; in fact the numerical 
estimate for the central charge, c ~ 1.60, is higher than the value c = 3/2, 
pertinent to a superposition of critical Ising and gaussian models [11]. Indeed 
the central charge seems to vary continuously from c ~ 1.5 near P to c ~ 2 
at T [10]. The system lacks conformal invariance [9], so one could consider 
suitable perturbations of the XY-Ising model as a starting point to study 
the vicinity of the point P [12]. Instead, recent numerical simulations on 
huge lattices [13] lead to two very close but separate transitions on the PT 
line. A possible solution could be [7] the addition of an antiferromagnetic 
coupling (J) term to the Coulomb gas Hamiltonian Hcg] for J ^ ^, the two 



transitions on the PT line separate with the KT one occurring at a lower 
temperature [7]. On the other side, by adding higher harmonics contributions 
to the potential [11, 14], the possibility of a merging critical point T will be 
provided here in the context of a twisted CFT approach [1], which extends 
the results of Ref. [11], so recovering the whole phase diagram [9, 10, 12]. 

2 m-reduction procedure 

In this Section we recall those aspects of the twisted model (TM) which are 
relevant for the FFXY model. We focus in particular on the m-reduction 
procedure for the special m = 2 case [1], since we are interested in a system 
with f/(l) ® ^2 symmetry. Such a theory describes well a system consisting 
of two parallel layers of 2D electron gas in a strong perpendicular magnetic 
field, with filling factor z/(") = ^^—^ for each of the two a = 1, 2 layers [1]. 
Regarding the integer p, characterizing the flux attached to the particles, we 
choose the "bosonic" value p = 0, since it enables us to describe the highly 
correlated system of vortices with flux quanta ^. Let us start from the 
"filling" u = ^, described by a CFT with c = 1 in terms of a scalar chiral field 
Q{z) = q — ip Inz + i X]n=^o ^'^~"' compactified on a circle with radius R^ = 
l/u = 2; here a„, q and p satisfy the commutation relations [a„, a„/] = n6ny 
and [q,p] = i. From such a CFT (mother theory), using the m-reduction 
procedure, which consists in considering the subalgebra generated only by 
the modes in Q{z) which are a multiple of an integer m, we get a c = m 
orbifold CFT (the TM), which is symmetric under a discrete Zm group and, 
for m = 2, will be shown to describe the whole phase diagram of the FFXY 
model. Its primary fields content, for the special m = 2 case, can be expressed 
in terms of two scalar fields given by: 

X{z) = \ {Q{z) + Q{-z)) , 0(^) = i {Q{z) - Qi-z)) ; (3) 

X{z) is Z2-invariant and describes the electrically "charged" sector of the 
new theory, while (f){z) satisfies the twisted boundary conditions (f){e^'^z) = 
—(piz) and describes the "neutral" sector [1]. The TM primary fields are 
composite vertex operators V (z) = Ux (z) ip (z) , where Ux {z) is the vertex 
describing its "charge" content and ip (z) describing the "neutral" one. In the 
neutral sector it is useful to introduce the two chiral operators ip (z) (ip (z)^ = 



_i ^. ^ta4>{z) . _j_ . ^ta4>{ z) .^^ -vi^ith only the first one obeying tlie boundary 
conditions. In a fermionized theory they correspond to two c = 1/2 Majorana 
fermions with Ramond and Neveu-Schwartz boundary conditions [1] and, in 
the TM, they appear to be not completely equivalent. In fact the whole 
TM theory decomposes into a tensor product of two CFTs, a Z2 invariant 
one with c = 3/2 and symmetry f/(l) ® Z2 and the remaining c = 1/2 
one realized by a Majorana fermion in the twisted sector. Furthermore the 
energy-momentum tensor of the Ramond part of the neutral sector develops 
a cosine term, T^ (z) = —j (90) — j^ cos (2v^0) , a signature of a tunneling 
phenomenon which selects out the new stable c = 3/2 vacuum; we identify 
such a theory with the one describing the FFXY model conjectured in Ref. 
[11]. 

3 FFXY phase diagram from TM model 

In this Section we will derive the FFXY phase diagram in terms of the RG 
flow which originates from perturbing our TM model with relevant operators. 
We observe that the limit /i — > in the Hamiltonian of Eq. (2) gives rise, in 
the continuum, to a CFT with two scalar boson fields f^*^*-* and with central 
charge c = 2. Now a good candidate to describe the FFXY model at criti- 
cality around the point T of the phase diagram is a CFT, with c = 2, which 
accounts for the full spectrum of excitations of the model: vortices, domain 
walls, and kinks. The role of the boundary conditions in the description of 
the excitation spectrum is crucial. In fact, by imposing the coincidence be- 
tween opposite sides of the square lattice, we obtain a closed geometry, which 
is the discretized analogue of a torus and gives rise, for the ground state, to 
two topologically inequivalent configurations, one for even and the other one 
for odd number of plaquettes. So the ground state on the square lattice maps 
into the ground state for the even case while it generates two straight domain 
walls along the two cycles of the torus for the odd case. Such a behaviour 
has to be taken into account by non trivial boundary conditions on the field 
(/9*^*^ at the borders of the finite lattice. To this aim, let (— L/2,0), (L/2,0), 
(L/2,L), {—L/2,L) be the corners of the square lattice C and assume that 
the fields y?'-*-' satisfy the boundary conditions cp^^^r) = ip^'^\r) for r G £nx, 
X being the x axis. That allows us to consider the two fields (p^^^^ and (p^"^^ on 
the square lattice £ as the folding of a single field Q, defined on the lattice 



£o with comers (-L/2, -L), (L/2, -L), (L/2, L), (-L/2, L). We can imple- 
ment a discrete version of the 2-reduction procedure by defining the fields 
X{r) = i(Q(r) + Q(-r)), $(r) = i(Q(r) - Q(-r)) (r G A), which are 
symmetric and antisymmetric under the action of the group Z2. The Hamil- 
tonian (2) can be rewritten in terms of these fields and, for h = 0, it gives 
rise, in the continuum, to the TM action -A. = J [^ (d^) + | (i9$) ] (fx. Let 
us now show how the phase diagram of the FFXY model can be described 
by the action: 



A 



- {dXf + - (9$)^ + /i cos (/5$) + A cos ( ^$ + 6 



(fx, 

(4) 



which embodies the higher harmonic potential term conjectured in Refs. [11, 
14, 7]. We assume the constraints 13"^ < Stt, which characterizes both the co- 
sine terms as relevant perturbations, and \5\ < tt/2 [15]. Thus the "neutral" 
sector is a two-frequency sine-Gordon theory that can be viewed as a defor- 
mation of a pure sine-Gordon one with the perturbing term A cos (/9$/2 -|- 5). 
The ultraviolet (UV) fixed point /i = 0, A = of the action (4) corresponds to 
the TM model with central charge c = 2, describing the fixed point T in the 
FFXY phase diagram. In order to study the RG fiow in the "neutral" sec- 
tor let us define the dimensionless variable r] = A/i~^ '^-(p/ ) )/[ '^-p j- when 
7] = Q the "neutral" sector reduces to a sine-Gordon model with a particle 
spectrum built of solitons and antisolitons and, for (3"^ < An, some breathers. 
Switching on the perturbation a confinement of solitons into states with zero 
topological charge takes place and packets formed by 2 of the original soli- 
tons survive as stable excitations for generic values of |5| < 7r/2. In the limit 
T) ^ 00 the 2-soliton evolves into the 1-soliton of the pure sine-Gordon model 
with n = 0. An unbinding phenomenon takes place in the case 6 = ±7r/2 for 
finite f] and the 2-soliton decomposes into a sequence of two kinks Ki. So 
the existence of an intermediate critical value r] = rjc'is required at which a 
phase transition takes place and the RG fiow ends into the infrared (IR) fixed 
point described by a CFT with central charge c = 1/2, the Ising model. The 
central charge of the full model (4) so changes from c = 2 of the UV fixed 
point to c = 3/2 of the IR fixed point, i.e. we recover early known Monte 
Carlo results [10]. Such an IR fixed point coincides with the U{1) ® ^2 sym- 
metric component of the TM model, which results then to properly describe 
[11] the fixed point P in the phase diagram of the FFXY model. 
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